% 生成新数据
clear,clc
diary on
s = RandStream('mt19937ar','Seed',235);
RandStream.setGlobalStream(s);
% original data
oriData = [459  362  624  542  509  584  433  748  815  505  612  452  434  982  640  742  595  706  593  680  926  653  164  487  734  608  428  1153  593  844  527  552  513  781  474  388  824  538  862  659  775  859  755  649  697  515  628  954  771  609  402  960  885  610  292  837  473  677  358  638  699  634  555  570  84  416  606  1062  484  120  447  654  564  339  280  246  687  539  790  581  621  724  531  512  577  496  468  499  544  645  764  558  378  765  666  763  217  715  310  851];
disp('原始数据的均值和标准差：')
mean(oriData) 
std(oriData)
fitPa = wblfit(oriData); % fit Weibull para
disp('原始数据的Weibull参数：')
disp(fitPa)

a = 600;
b = 4;
fprintf('a = %i, b = %i, Weibull 随机数：\n',int32(a),int32(b))
out = round(wblrnd(a,b,15,10));
disp(out)
disp('新生成数据的均值和标准差：')
mean(out(:))
std(out(:))
diary('data.txt')
diary off